DESY 12-042 
ZMP-HH/12-4 

Anisotropic hydrodynamics, holography and the chiral magnetic effect 



(N 

o 

(N 

1—5 

o 



Ph. 
CD 

(N 
> 

in 

(N 

en 

o 



X 



Ilmar Gahramanov, Tigran Kalaydzhyan, and Ingo Kirsch 

DESY Hamburg, Theory Group, Notkestrasse 85, D-22607 Hamburg, Germany and 
Zentrum fur Mathematische Physik, Universitdt Hamburg, Bundesstrasse 55, D-20146 Hamburg 

(Dated: February 27, 2013) 

We discuss a possible dependence of the chiral magnetic effect (CME) on the elliptic flow coeffi- 
cient V2- We first study this in a hydrodynamic model for a static anisotropic plasma with multiple 
anomalous f/(l) currents. In the case of two charges, one axial and one vector, the CME formally 
appears as a first-order transport coefflcient in the vector current. We compute this transport coef- 
ficient and show its dependence on V2. We also determine the CME coefficient from first-order cor- 
rections to the dual anti-de Sitter background using the fluid-gravity duality. For small anisotropies, 
we find numerical agreement with the hydrodynamic result. 

PACS numbers: 11.15.-q, 47.75. -ff, 11.25.Tq, 12.38.Mh 



I. INTRODUCTION 



In the last couple of years the chiral magnetic ef- 
fect (CME) has attracted much attention as a candi- 
date for the explanation of an experimentally observed 
charge asymmetry in heavy-ion collisions, as seen by the 
STAR [H, PHENIX and ALICE collaborations. 
The CME is a hypothetical phenomenon which states 
that, in the presence of a magnetic field B, an electric 
current is generated along B in the background of topo- 
logically nontrivial gluon fields [j, Q . Analogous effects 
were found earlier in neutrino [6|, electroweak [7| 
condensed matter physics . Lattice QCD results 
suggest the existence of the effect, although the magni- 
tude of the CME-induced charge asymmetry may be too 
small to explain the observed charge asymmetry jl^ . 

In a recent experiment, the charge separation is mea- 
sured as a function of the elliptic flow coefficient V2 ^^ ■ 
The data is taken from (rare) Au+Au collisions with 
20 — 40% centrahty but different V2- In this way V2 is 
varied while at the same time the number of participat- 
ing nucleons (and therefore the magnetic field) is kept al- 
most constant. The plots in [T^ suggest that the charge 
separation is proportional to V2- If this holds true, the 
charge separation will depend on the event anisotropy. 

In this paper we address the question of whether and 
how the CME depends on the elliptic flow V2- We study 
this both in hydrodynamics and in terms of a holographic 
gravity dual. The hydrodynamical approach to the CME 
and CME-related phenomena was proposed in p^| - |20| . 
There, the CME appears in form of a nonvanishing trans- 
port coefficient in the electric current, j = hbB, which 
measures the response of the system to an external mag- 
netic field [ij, [21| . In [23l , the chiral magnetic conduc- 
tivity in an isotropic fluid was determined as 



KB = Cn5 1 



+ P 



(1) 



The first term is the standard term for the CME and 
depends only on the axial anomaly coefficient C and the 
axial chemical potential /X5. The second term propor- 
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FIG. 1. Sketch of the time evolution of the momentum 
anisotropy Ep (based on [22l). The small figure shows the 
orientation of Pl and Ft with respect to the reaction plane. 



tional to the factor ^^ depends on the dynamics of the 
fluid and has a chance to depend on V2 in the anisotropic 
case. 

In the flrst part of the paper we study this in a hydro- 
dynamic model for an anisotropic fluid with multiple 
anomalous U{1) charges (This model extends those in 
J23H25I). We compute the CME coefficient kb and ex- 
press the result in terms of the momentum anisotropy Sp 
1261 defined as 



{Pt - Pl) 
{Pt + Pl) 



(2) 



where Pt and Pl are the pressures in the plane transverse 
to the beam line (In our conventions the indices L and T 
refer to the longitudinal and transverse direction with re- 
spect to an anisotropy vector w^ normal to the reaction 
plane, see Fig. [IJ. A sketch of Sp as a function of the 
proper time r is shown in Fig. [1] Ep describes the build- 
up of the elliptic flow in off-central collisions. Our model 
describes a state after thermalization with unequal pres- 
sures Pt i^ Pl- At freeze-out Sp roughly equals V2^ and 
we flnd that for small anisotropies the CME-coefflcient 
Kb increases linearly with V2- 



In the second part of the paper we perform a holo- 
graphic computation of kb hi the dual gravit y n iodcl. A 
similar com put ation was previously done in [20[ for the 
STU model [231 j ^ string-theory- inspired prototype of an 
(isotropic) anti-de Sitter (AdS) black hole solution with 
three C/(l) charges. Other holographic approaches to the 
CME can be found in (H-IH. 

In the anisotropic case, we first need to construct an 
appropriate gravity background. As an ansatz, we choose 
a multiply charged AdS black hole solution with some ad- 
ditional functions w^ and wt inserted which will make 
the background anisotropic and £p-dcpcndent. Since ana- 
lytical solutions for charged anisotropic backgrounds are 
notoriously difficult to find, we will use shooting tech- 
niques to find a numerical solution. Other AdS back- 
grounds dual to anisotropic fluids are constructed in |39l - 

Eii. 

As the AdS solution in [3^, the background is static 
and does not describe the process of isotropization. Even 
though such models have some limitations [39| , they are 
nevertheless useful for the computation of transport co- 
efficients. We show this, following [201, by determining 
Kb from the first-order corrections to this background us- 
ing the fluid-gravity duality ^^. For small anisotropics, 
we find numerical agreement with the hydrodynamic re- 
sult for Kb- Other (dissipative) transport coefficients 
in strongly-coupled anisotropic plasmas are discussed in 

The paper is organized as follows. In Sec. [TTl we re- 
view the hydrodynamics of an anisotropic relativistic 
fluid with several U{1) charges and triangle anomalies. 
We then compute the vortical and magnetic conductivi- 
ties of such a fluid by extending the method of Son and 
Surowka [IJI to the anisotropic case. In Scc. lIIII wc con- 
struct the dual gravity background and present a numer- 
ical solution for its gauge field and metric functions. In 
Sec.|lV]we use this background to perform a holographic 
computation of the vortical and magnetic conductivities. 



In anisotropic relativistic fluids, the hydrodynamic 
equations are again given by ([3]) but the stress-energy 
tensor T^'^ and U{1) currents j°^ now have the more 
general forrrl^ 



Tt"" = (e + PT)u^'u'' + Prg^'" - Av^'v'' + t^' 

■afji 



(4) 

J-" = p°M^' + y'^t' , (5) 

where e is the energy density, p" are the f/(l) charge 
densities, A = Pt — Pl, and Pt and P^ denote the 
transverse and longitudinal pressures, respectively |23l - 
[25|. 5^1/ is the metric with signature (—,+,+,+). t'^'^ 
and 1'°''^ denote higher-gradient corrections, for which we 
require u^t^'^ = and u^ji/"'' = 0. 

The four-vectors u^^ and v^^ describe the flow of the 
fluid and the direction of the longitudinal axis, respec- 
tively. The vector v^ is spacelike and orthogonal to u^ , 

w^M^ = -l, w^w'' = l, u,.«^=0. (6) 

It is convenient to define the proper time t by S"^ In r = 
v^'d^v" 11]. In the rest frame of the fluid, u^' = (1, 0, 0, 0) 
and v^ = (0,0,0,1), the stress-energy tensor becomes 
diagonal, 

e 

Pt 

Pt 

Pl 

In conformal fluids, the stress-energy tensor is traceless, 
T^p = 0, and e = 2Pt -I- Pl- Clearly, the isotropic case 
corresponds to equal pressures Pt and Pl, P ~ Pt ~ Pl ■ 

For simplicity, we restrict to the case of a single charge 
in Sees. HIA] and EbI n = 1. In Sees. HC] and ED] we 
generalize our findings to arbitrary n and discuss the case 
n = 2, which is relevant for the CME. 



Thermodynamics of an anisotropic fluid w^ith 
chemical potential (n = 1) 



rpfj,!. 



(7) 



II. HYDRODYNAMICS OF ANISOTROPIC 
FLUIDS WITH TRIANGLE ANOMALIES 

The hydrodynamic regime of isotropic relativistic flu- 
ids with triangle anomalies has been studied in [14| - [l9| . 
and much can be taken over to the anisotropic case. Such 
fluids typically contain n anomalous U{1) charges which 
commute with each other. The anomaly coefficients are 
given by a totally symmetric rank-3 tensor C°^'^. The 
hydrodynamic equations are 



d^T'^" 



pau\-a^ 



^^d'"' = C^'^E" ■ B 



(3) 



where E"!" = P^^^^w^,, P''^ = ie^'^^^u^P^^ (a = 1, ..., n) 
are electric and magnetic fields, and P° = d^A^j — 9,yA" 
denotes the gauge field strengths. As in [ij], we expand 
the constitutive equations for T^"^ and j^ up to first or- 
der, taking A'^ ^ 0{p^) and P^^ ^ 0{p). The gauge 
fields A° are nondynamical. 



Hydrodynamic models for an anisotropic fluid (without 
chemical potential) have been studied in [23l - |25| . Follow- 
ing these works, we derive some thermodynamic identi- 
ties, now for the case of a fluid with a chemical poten- 
tial pi. 

These identities can be found by computing the quan- 
tity /o = Uudfj.T^^" + fid^j^ at zeroth order. Since the 
right-hand side of (j3]) can be dropped at order 0{p'^), we 
have /o = 0. Using d^{su'^) = 0, we get 

u^d^.Ti"' = -uf'df.e - (e + PT)d^,u^' - Au^d" Inr 



A 



S T 



fidf^r = M(5^P)"^ - ^W'd^s . 



(8) 
(9) 



^ The symmetries allow in principle for more general currents 
ja/i _ pa^^fj, _|_ (.tij^M _|_ j^ap ^j^^j^ some Coefficients c°. Here we 
switch off all the 'electric' background currents, c" = 0. 



As in [241 ■ we consider a generalized energy density 
e = e(s,/9, r), which depends not only on the entropy 
density s and particle density p but also on the new vari- 
able T. Its differential is 



de 



with 
ds 



ds 



dp 



dp 



dr, (10) 



.r"^' i^pJs^r'- 



A 



(11) 



The temperature and the chemical potential are defined 
in the usual way. If we also impose [de/dr)^ = — A/r 
and substitute PU]) into (|S]), then /q = implies the 
following thermodynamical identities for an anisotropic 
fluid: 



e + Pr- 


= Ts + pp , 


(12) 


dPr -- 


= — dr + sdT + pdp , 

T 


(13) 


de-- 


A 
= Tds + /irfp dr , 


(14) 



in agreement with [2J| for /.j = 0. 



B. Vortical and magnetic coefficients (n = 1) 



We now discuss corrections to the L'^(l) current 
jfj- ^ jiM (^ = i^^ Iyy anisotropic fluids the transport co- 
efficients are usually promoted to tensors such that one 
should consider first-derivative corrections of the type 



^^^ = (e.)^C.''-f(eB)%S^ 



(15) 



where w^ = ^e^P^^u^dpU^ is the vorticity, and B^ is 
an external magnetic field. In Landau frame u^v^ = 
and therefore u^{S,ojYu^'^ = (and similar for {£,bYu)- 
This is satisfied e.g. for (^w)'';y = ^^^^i/, since u^uj^ = 
(We do not consider other components of ^^^ here). We 
therefore restrict to consider corrections of the type 



i^Lu^ + £,bB^^ ., 



(16) 



as in the isotropic case [lj|. Our goal is to compute the 
vortical and magnetic conductivities ^(^ and S,b- These 
transport coefficients can be found by assuming the exis- 
tence of an entropy current s^ with a non- negative deriva- 
tive, S^s'' > 0. The computation closely follows that 
of [i|. 

The hydrodynamic Eqs. (j31) imply that the quantity 



h = u,^^T^''' + pd^f + E^'v^ - pCE^'B, 



(17) 



vanishes at first order, /i = 0. Substituting the explicit 
expressions for the stress-energy tensor and U{\) currents 



into /i and using the thermodynamical identities (J12p 
and (fT4l). we find 



d^,[s 



— I 
T 






-C^E-B, 
T 



(18) 



which is exactly the same equation for the entropy pro- 
duction as in the isotropic case [ij]. 

In the following, we will need the identities 



d^Lj'' = 



P 



-uj^id^Pr - Ad^ In T-pE^), (19) 



a^B^ = -2lu^'E^ - 



TPi 



-{O^^Pt - Adf^lnr - pE^) , 



which we derived from ideal hydrodynamics in Ap- 
pendix [X] In deriving these identities we assumed that 
the fluid satisflcs 



d^v'' = , 



'c'„A = 0. 



(20) 



The first equation is basically a "continuity equation" 
for the vector v^^. There are no sources for the genera- 
tion of anisotropy. The second equation imposes an or- 
thogonality relation between the gradient of the pressure 
difference A = Pt — Pl and v^. 

As in |14| , we assume a generalized entropy current of 
the form 



— I 

T 



Duji" + Db^^ 



(21) 



where S,^, ^b, D, and Db are functions oi T, p and r. 
We now compute S^s'*, using (fT8|) and ([T9|) and impose 
dfj.s^ > 0. Since the coefficients in front of w^, i?^, tu^E^^ 
and Efj^B'^'' inside d^s^ can have either sign, we require 
them to vanish and obtain the following four differential 
equations: 



d^D 



-^^ {d^Pr - Ad^ In r) - ^^8,, ^ = , (22) 
- -^{d^PT ^ Ad.lnr) ^ ^Bd,^ = , (23) 



e + Pt T 

e + Pr T T 



(24) 
(25) 



For A = 0, these equations reduce to those in the 
isotropic case |14| . 

In Appendix M we solve (I121)-(Iin|) for D, Db, C^ 
and S,B- As a result, we find the vortical and magnetic 
conductivities 



^u^^Cip'-- 



,2 2 pp^ 



3e + PT 



0{T^), 



(b=C[p-\-P^]+0{T^), (26) 



where 0{T^) denotes terms proportional to T^, see (|B14[) 
in Appendix [BJ These terms are related to gravitational 
triangle anomalies [I^, |4g and may, in the anisotropic 
case, depend on the proper time r. In the absence of 
gravitational anomalies, which we do not discuss in this 
paper, the conductivities do not depend on r. Apart 
from these changes in 0{T^)^ the relations have the same 
form as in the isotropic case but with P replaced by the 
transverse pressure Pt- 



C. Multiple charge case (n arbitrary) 

The generalization of the previous computation to a 
fluid with multiple anomalous U{\) charges is straightfor- 
ward, and we only state the result here. The corrections 
v°-^^ of the currents j"^ in ^ are 



(AAA) 



with [terms of order O(T^) ignored] 



bc.,b,,c 2 ar^bcd^^ f^'^ f^ 



3 e + Pt 



(:ab 



riabc,,c ^ „a/-ibcd /^ ^^ 



(27) 

(28) 
(29) 



These are simple generalizations of the corresponding 
conductivities in the isotropic case [ij, ^^ . 



D. Chiral magnetic and vortical effect (n = 2) 

Physically, the most interes ting case is that involving 
two charges (n = 2) [la, [l3, 120] ■ The chiral magnetic 
effect [j| can be described by one axial and one vector 
U{1), denoted hy U (!) a ^ U {l)v ■ A convenient notation 
for the gauge fields and currents is (a, b, ... = 1, 2) 



A"^ — 41 



A^ - A^ 



(30) 



Let us now derive the chiral magnetic and vortical 
effects from ([28|) and ([29| . C— parity allows for two 
anomalous triangle diagrams, (AAA) and (AVV), shown 
in Fig. m while diagrams of the type (VVV) and (VAA) 
vanish. Accordingly, the anomaly coefficients are 

^121^^211^^112^0^ {VAA) 

6*222 ^ Q ^ (yvv) 

C"VO, {AAA) 

^122 ^^221^ ^212 ^0_ ^J^yy^ 



(31) 



The hydrodynamic Eqs. ([3]) then imply nonconserved 
vector and axial currents 

5^J^ = -i(C2i2F4,F^^- + C22iF^^,F^A.-) , 

<9pJ-r - -\{C^^'F^^F^^^ + C^-'^Fl^F^n , (32) 



(AVV) 





FIG. 2. Anomalous diagrams corresponding to C'^^^(left) and 
to C^^^ = C^^^ = C^^^ (right). Dashed (wavy) lines denote 
the axial (vector) currents/fields. 



afii^ 



where we rewrote £>> ■ B" = -jF^^^F"'"' (with F" 

Ip/iiypo-pa ^ 
2 ^ ^ pa )■ 

To restore conservation of the vector current, we add 
the (topological) Bardcen term to the boundary theory. 



Sb^cb I d'^xe^'^^PA^AlF, 



V 

\p- 



(33) 



Combining the corresponding Bardcen currents 

3''B^CBe^'''"\AXF^^~2AtFl), 
jlB^^'Be^^'^'AlFl, 
with the vector and axial currents, 

J —J^Jb' -lb ~ -JR ^ J5,B ^ 

we obtain the anomaly equations 



(34) 
(35) 

(36) 



/^iii 

Pi „-'M _ pA rpAap 



C^ 



CB I'.b *" , 



The electric current j''' is conserved if cb = —C^"^"^ /2. 
Setting C"! = 6*122 '^ ^^3^ ^j^g hydrodynamic Eqs. © 
become 



9^.?^^ = CE-B + {C/3)E5 ■ B5 . 
Using the derivative expansion 



(37) 



(38) 



where k^^ = ^2: '^s = Cs^ ^'^'^ ^^5.5 = ^s^, wc obtain from 
and (P5|) the conductivities (/is = /i^, 11 = fi^) 

P 



Kuj = 2C/.J5 ( /i - 
Kb = C*/i5 I 1 - 



e + Pr 
MP 



e + Pr r 



n , . ^ PP 



2e + PT 



9 M5 



Mi 



3/^2 



(39) 



There are analogous transport coefficients in the axial 
current j^ [20[. The axial fields E^^^ and Bs^ are not 
needed and can now be switched off. The first term in Kg 
and Kun i^B = C^5 and k,^ — 2Cun^, is the leading term 
in the chiral magnetic (CME) [J, |5| and chiral vortical 
effect [431 1 respectivelyo They are in agreement with 
those found in the isotropic case [la, [l3, 1201 ■ The second 
term proportional to p/(e + Pt) actually depends on the 
dynamics of the fluicO and therefore on Sp . 

The dependence of kb on Ep can be made more visible 
by introducing an average pressure P ~ {2Pt + Pl)/3 
such that e = 3P. Assuming Sp to be small (see Fig. [TJ, 
we expand the CME-coefficient kb to linear order in Sp, 



KB « C/i5 1 - 



fip 



P 



-? 



(40) 



At freeze-out the elliptic fiow coefficient V2 ~ £p/2 j26j . 
For small momentum anisotropics, the CME thus in- 
creases linearly in f 2 . 



equations. 

Gain — (^gMN = Tj 



^mF'^^'p = - 



MN 



(44) 



be PMNKL Tj^b 



87=5 
where the energy-momentum tensor Tmn is 



FIinF'kl , (45) 



Tmn = -2 ( FI^rF-^n + -,9mnF^srF''^'' ) ■ (46) 



A. AdS black hole with multiple U(l) charges 

A gravity dual to an isotropic fluid (e = 3P) with mul- 
tiple chemical potentials p,a (a = l,...,ri,) at finite tem- 
perature T is given by an AdS black hole solution with 
mass TO and multiple f/(l) charges q"'. In Eddington- 
Finkclstein coordinates, the metric and U{\) gauge fields 
of this solution are 



III. FLUID-GRAVITY MODEL 

In this section we construct the gravity dual of a static 
anisotropic plasma with diagonal stress-energy momen- 
tum T^i, ~ diag(e, Pt, Pt, Pl) and charge densities p". 

We start from a five-dimensional t/(l)" Einstein- 
Maxwell theory in an asymptotic AdS space. The action 
is 



1 



167rG5 



(Px^ 



2A 



T^a jpaMN 



(41) 



Oabc ^PKLMN Aci jpb jpc 
a I — -^ ^P^ KL^ MN 



where A ~ —6 is the cosmological constant. As usual, 
the C/(l) field strengths are defined by 



^MN = dn^N ~ QnAIj , 



(42) 



where M,N,... = 0, ...,4 and a = l,...,n. The Chern- 
Simons term A A F A F encodes the information of the 
triangle anomalies in the field theory [ij]. In fact, the 
Chern-Simons coefficients Sabc ai'c related to the anomaly 
coefficients Cabc by 



C'abc = Sabc/i^T^G5) 



(43) 



The corresponding equations of motion are given by 
the combined system of Einstein-Maxwell and Maxwell 



ds^ = -f{r)dt^ + 2drdt 
A"- = -Al{r)dt, 



r^dx^ . 



(47) 



where 



/(o = ^^-^ + E 

a 



m' 



^S(0 = a4 + 



2r2 



(48) 



The constants ^J^ can be fixed such that the gauge 
fields vanish at the horizon. In case of a single charge 
{n= 1), the background reduces to an ordinary Reissner- 
Nordstr0m black hole solution in AdS^ [4^. 

The temperature T and chemical potentials /^° of the 
fluid are defined by 



T 



2^ 



/'(r+)_2r6-^Jg,) 



47r 27rr^ 



(49) 
(50) 



where r+ is the outer horizon defined by the maximal 
solution of f{r) = 0, and r^o indicates the location 
of the boundary. The temperature of the fluid is the 
Hawking temperature of the black hole and is computed 
from the surface gravity k = V^^aTIxK^^xI I r+ 1 where 
Ixl = (—X^Xm)^^^'^'' is the norm of the timelike Killing 
vector x^' = C [here |x| = //Wl- 



^ 1^5,3 represents another effect, which we added for completeness, 
but it seems not to be reahzed in heavy-ion collisions. 
In |15l | this term was considered as a one-loop correction in an ef- 
fective theory and ie + P)/ p was interpreted as the corresponding 
infrared cutoff in the energy/momentum integration. 



B. Anisotropic AdS geometry with multiple U(l) 
charges 

We now construct a solution for an anisotropic fluid 
(e = 2Pt + Pl)- An ansatz for an anisotropic AdS black 



hole solution is given by 

ds^ ^ -f{r)dt^ + 2drdt 

+ r {wT{r)dx +WT{r)dy + WL{r)dz ), 
A" ^-A^{r)dt. (51) 

The anisotropics arc realized via wt{t) and WL{r), which 
are functions of the momentum anisotropy Sp as defined 
in ©, 



(Pt 



Pt + Pl 



(52) 



In the isotropic case {Sp = 0), these functions are required 
to be one, wrir) = WL^r) = 1, and the background re- 
duces to the AdS black hole geometry (|T7l) . 

An analytical solution of the type ([51]) is difficult to 
find, and we resort to numerics in the next subsection. 
For this, we need to know the solution close to the bound- 
ary. An asymptotic solution (?' — > cx)) is given by the four 
functions 



2r2 



fir)y 



m 
r 



i+E^+^(^-^)^ 



w. 



WT{r) = 1 + ^ + Oir-"") 



(4) 
T 

„(4) 



WLir)^l + ^+Oir-^)., 



(53) 



where w 



(4) 



-2w, 



(4) 



-mC,/2, /ij^ = const., and C, is 



related to the momentum anisotropy £„ by 



c = 



2e„ 



Ep -I- 3 



(54) 



The functions wt{t) and WLir) have been intro- 
duced in view of the structure of the anisotropic fluid 
stress-energy tensor. More precisely, in ([53| we fixed 



the 



,(4) 



coefficients wL and w\ ' such that the fluid 



(4) 



stress-energy tensor is of the diagonal form ([7]), T'^'^ = 
diag(e,PT, -Pt,-Pl) with e = 2Pt + Pl- Computing the 
stress-energy tensor in the standard way from the asymp- 
totic solution (j53p via the extrinsic curvature, see e.g. 
[49| , we find the transverse and longitudinal pressures 



Pt = 



Pt, 



m — Aw, 



(4) 



4w 



(4) 



167rG5 

„(4) 



m(l + C) 
IGTrGs 



m - ^w)^' m(l - 2C) 



167rG5 



167rG5 



(55) 
(56) 



Note that if ([5^ holds true, the pressures Pt and Pl 
satisfy ((52|) . Likewise, the charge densities are 



167rG5 



(57) 



From these relations, we find the useful identity 

^^ ^ V3g° 

e + Pr 4m(l + iC)' ^ ' 

which we will need later. 

Numerical solution 

We now use shooting techniques to solve the system 
of ordinary differential equations (ODE) which follows 
from the equations of motion pi)) and (PSj) upon substi- 
tuting the ansatz ([5T|) . The idea is to vary the metric 
and gauge fields at some minimal value r+ in the radial 
direction, integrate outwards and find solutions with the 
correct asymptotic behavior (j53|) . A similar method was 
previously applied in [4l|. 

We first need to study the asymptotic solution near r_|_ 
and near the boundary at roo 3> r+ (we choose Too = 50 
in our numerics). We define r_|_ by the maximal solution 
of 



/(' + ) = 



(59) 



and use scale invariance to set r+ = 1. We then ex- 
pand the functions in the metric and gauge fields near 

r+ in powers of the parameter e = — 1^1 and 

substitute them into the equations of motion. In this 
way, we find that the only independent variables are 
{f'{r+),WT{'i'+),WL{r^),w'j^(rj^)} since the gauge field 
parameters A^irj^^) can be set to zero using gauge in- 
variance, Ag(r+) = 0. The other parameters at r+ 
can be expressed in terms of these four parameters, e.g. 
WT{r+) = WT{r+)w'L{r+)/wL{r+). 

The near-boundary solution is given by (j51|) with (|53p 
and is parameterized by the values {C,,m,q°' , ^'^). The 
final set of data is summarized in the following table: 



r = r+ = 1 


r = rao:>r+ 


Ag(r+)=0 


MSo 


/('^+) = 


Kroo) 


/'(r+) - fixed 


Al'{r^) 


WL{r+) = var 


WL{roc) 


wt (?"+ ) = var 


WT{r oo) 


w'L{r+) =var 





Parameters not listed are related to those in the table by 
the equations of motion. 

To integrate the equations we proceed as follows. We 
fix C and vary three parameters at 7-+, namely wt (''+), 
WL{r+) and i(j^(r+), by choosing a grid with suitable 
number of sites (in our case 20'^ — 40"^). The value f'{r+) 
can be thought of as the temperature of the system and 
will simply be fixed to some value. It turns out that the 
form of the functions WL,T{r) does not depend on this 
parameter. For each site in the grid we numerically solve 
the system of ODEs and determine the pair (m, g") from 



the known asymptotics of ^o'(r 



d) and /(r 



.)■ 



This ensures that the analytical and numerical values for 
these quantities coincide. 

We then calculate the combined residual 



reso 



[wT{r+),WL{r+),w'^{r+))] 
{w*{r^) - wlir^))^ + iw*ir^) - w*T{r^))^ (60) 



„# 



where wf^ r^{r oa) are the numerical values, and w^ t(''oo) 
are the analytical values given by ([53]) . We interpolate 
the residual by a piecewise linear function and find its 
global minimum by the simulated annealing method [50[ . 
The result of the minimization is shown in Fig. [31 which 
depicts numerical plots of /(r), Ao(r), wrir) and WL(r) 
for n = 1. 

-Ao[r],f[r]/lo3 
3.5 




FIG. 3. Numerical plots of /(r), Aoir), wrir) and WL{r) for 
^ = 10 (r+ = 1). We get WLir+) = 12.42. 

We conclude this section with a comment on r+. In 
the isotropic case, r_|_ is simply the size of the hori- 
zon of the AdS black hole geometry. For nonvanishing 
anisotropics and vanishin g U (V) charges, a naked sin- 
gularity was found at r+ [39[, implying that the static 
background does not exist indefinitely. The singularity is 
mild in the sense that there is a notion of ingoing bound- 
ary conditions and possible instabilities are absent at the 
linear level in the anisotropy parameter [39|. This be- 
havior may persist even for nonvanishing C/(l) charges, 
even though it was difficult to see the singularity in our 



numerics, c/. Figure SI Despite this subtlety, we show 
in the next section that, at least for small anisotropics 
where the bulk geometry approximates a black hole so- 
lution, the singular geometry may be used to compute 
some transport coefficients of the fluid. 
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20 



FIG. 4. Numerical plots of {RMNpqf for C = 10, g / (red), 
C = 10, g = (orange), and C, — 0, q — (blue). 



IV. HOLOGRAPHIC VORTICAL AND 
MAGNETIC CONDUCTIVITIES 

We will now compute the chiral vortical and magnetic 
conductivities ^^ and ^^^ from first-order corrections to 
the numerical AdS geometry (j5ip using the fluid-gravity 
correspondence j43 |. 



A. First-order corrected background 

In order to become a dual to a multiply charged fluid, 
the AdS geometry ([5T|) must be boosted along the four- 
velocity of the fluid w^ (/.t = 0, ..., 3). The boosted version 
of (ISTI) is 



ds^ = {r'^WT{i')Pfj,i, — .f{r)u^Ui,^ dx^dx'^ — 2u^dx^dr 
— r'^{wT{r) — w L{r))v ^j,v,jdx^ dx^ , 



A'^ = {A'^^{r)u,,+Al)dx 



(61) 



where Pf"" == g'"'' + u^'u'' , and /(r), Ag(r), wrir) and 
WLir) are numerically known functions. As in hydro- 
dynamics, the four- vector v'^ determines the direction 
of the longitudinal axis, cf. Sec. 2. Following [ij, [20{ . 
we have formally introduced constant background gauge 
fields Afi^ to model external electromagnetic fields, such 
as the magnetic fields B"'^ needed for the chiral magnetic 
effect. 

The transport coefficients ^° and S,]^ can now be com- 
puted using standard fluid-gravity techniques [421. We 
closely follow [ij, [20, [51| , in which these transport coef- 
ficients were determined for an isotropic fluid with one 



and three charges {n = 1,3). We work in the static frame with 
u^ = (—1, 0, 0, 0), w^ = (0, 0, 0, 1), and consider vanishing 
background fields A°^ (at x^^ =0). The transport coeffi- 
cients ^" and ^^'' measure the response of the system to 
rotation and the perturbation by an external magnetic 
field. We therefore slowly vary the velocity u^ and the 
background fields Af^ up to first order as 



u^ = {-l,x-'d,u,), Al^iO^x^d^Af). 



(62) 



We may also vary m and q in this way, but it turns out 
that varying these parameters has no influence on the 
transport coefficients ^° s-^id £,'j^. 

Because of the dependence on x^^ , the background (pT|) 
is no longer an exact solution of the equations of mo- 
tion. Instead with varying parameters the solution (j6ip 
receives higher-order corrections, which are in this case 
of first order in the derivatives. 

An ansatz for the first-order corrected metric and 
gauge fields is given by 



ds^ = (-/(r) + gtt) dt^ +2(1+ gtr) dtdr 

+ r'^{wT{r)dx'^ + WT{r)dy^ + WL{r)dz'^) 

+ (jijdx^dx-' — 2x^ dijUidrdx'' 

+ 2 ((/(r) - r^) x^d.u, + ~gu) dtdx' , 

A'^=[-A^„{r)+At)dt 

+ (Ag(r)a;''a,u, + x^d.A^ + i^) dx' , 
where the first-order corrections are denoted by 

hiN = 9MN{r) , A^/ = Alj{r) . 
As in [5l| , we work in the gauge 



0, 9. 



rfi 



A° = 0. 



Y^ Oil = . 

1=1 



(63) 



(64) 



(65) 



The first-order corrections can be obtained by substi- 
tuting the ansatz (|63p into the equations of motion 
(|33|) and (|^5|) . The computation is straight- forward but 
lengthy and has been shifted to Appendix [C] [we set 
fi'^ = A^{roo) — there, see Sec. Ill CI for a discussion]. 
As a result, we find the following corrections: 



gtr = gtt = 
9tt{r) = f{r) 



At 



0, 



dr' 



AUr) 



r+ 

dr' 



WL{r'y/'^r'{f{r')f 
dr"I{r")-WLir+y/'r+rir+)C, 
1 



(66) 



o r'f[r')wL{r'YI^ 
C,r+A^'{r+)wL{r+f'^ + r'~gu{r')A'^^' {r') 



/(r)=^4<(r)(Q^(r)-g^(r+) 

a=l 



C 



c{r+) 



4c(r+) 
X (lSa,cA^oir+)AUr+)A^oir+)e'"' {d.u,) 



+ -SabcA^,{r+)Al{r+)^=^ {d,Al) y 



r+(/'(r+)-4Ea^S(''+K'(^+))' 



and r_|- as defined around (j59p [gij can be obtained by 
solving (|C2p in Appendix [C] but will not be needed here]. 



B. Holographic conductivities 

On the boundary of the asymptotic AdS space (|63p . 
the metric and gauge fields couple to the fluid stress- 
energy tensor and U{\) currents, respectively. Holo- 
graphic rcnormalization [52| provides relations between 
these currents and the near-boundary behavior of their 
dual bulk fields. For the magnetic and vortical effects, 
we need the U{1) currents .7'°^, which are related to the 
bulk gauge fields A°^ by [S^, l53| 



lim 



BttG. 



-V^^Alir). 



(67) 



Expanding the solution in - and substituting only the 
corrections A" we get the currents 



r^oo SttGs 

= Y^^V'"' {Q>+) + r+AS'(r+)a) . (68) 

Note that, in the isotropic case {wl = 1, Pt = Pl = 
P), the prefactor of the second term of (pS]) is simply 



r+<(r+)c(r+) = ^g" 



(69) 



as can be seen by substituting the Reissner-Nordstr0m 
solution (P5| into the left-hand-side of this equation. In 
the anisotropic case, we need to show that 



r+<(r+)c(r+).«;i(r+)-i/2 



1 



\/3g° 

4m ' 1 + k 



(70) 
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FIG. 5. Values of WL{r+) as a function of the anisotropy (. 
The numerically determined values for WL{r+) lie on the solid 
curve, which represents the function (1 + \C)^ ■ 



which, by ([55)1 . is equivalent to p"-/{€ + Pt). This equa- 
tion holds in particular if the first and second factors 
on both sides agree individually. The first factors corre- 
spond to (|69| . vifhich is expected to hold, at least approx- 
imately for small anisotropies (. The second factors are 
identical if i(;i(r+,C) = (1 -|- |C)^- We find numerically 
(for n — 1) that WL{r+) indeed satisfies this equation, 
see Fig. [5] Thus ([70)1 holds numerically, at least in the 
limit of small C- 

Comparing (|68p with the general expansion 



= C ^e'^P^'^u^dpU, + ^Be'^P'^^^u^dpA", , (71) 
we finally obtain the coefficients 



C = TTTZTT I 5"''''mV' - \ ^° 5''"=VV'/u'' I , (72) 



^B 



IGTrGs 

4 
167rG5 



3e + PT 



gabc^c _ }__[^sbcd^c^d 



2e + PT 



(73) 



the action to obtain a conserved chiral charge [30[. This 
charge is however only gauge-invariant when integrated 
over all space in homogeneous configurations. 

In AdS black hole models of the CME, one usually 
introduces a chiral chemical potential dual to a gauge- 
invariant current, despite it being anomalous |20l . |3l|. 
The prize to pay is the appearance of a singular bulk 
gauge field at the horizon, a phenomenon which seems to 
be generic in AdS black hole models of the CME. 

Careful holographic renormalization shows that, in the 
presence of Chern-Simons terms, there is an additional 
term on the right-hand side of (p7| [53j . This term is of 
the form 






>-'a. 



be 



SttG. 



'A^°J{x)dpA(Pj{x), (74) 



where Aail{x) are the Oth-order coefficients in a - ex- 
pansion of the bulk gauge fields Aaij,{r,x). In ([65|) we 
expanded the background gauge fields ^° around zero 

and set A'^^ = ^.%^Uv = 0. This allowed us to ignore 
terms in ([57)) coming from ([71)1 (at least to first order in 
the derivatives). 

Problems arise if /i^ ^ 0. To sec this, let us restrict 
again to two charges {n = 2) as in Sec. Ill CI and define 
axial and vector gauge fields by A^ = Aj^ and A^ = A^ . 
Then j^ = J2 gives rise to additional contributions of the 
type 



i^De'-'"^A;^(«)(x)Fj/)(x), 



(75) 



which are forbidden by electromagnetic gauge invari- 
ance [301, unless A^^ '{x) = 0. However, in general 
A^^ '{x) = fif'Ui, (at a; = 0) with some constant nf. 
We should thus set jj."^ = [Note that this does not im- 
ply fi5 = A^(roo) — A^{r+) = 0]. This corresponds to a 
nonvanishing gauge field at the horizon, as noticed also 
in [13,1311. 



V. CONCLUSIONS 



with fi°- EE ^o(r+) [since ^g(oo) = 0]. Using the re- 
lation (|43|. we find that the holographically computed 
transport coefficients (|72p and (|73p coincide exactly with 
those found in hydrodynamics, (|28p and 



C. Subtleties in holographic descriptions of the 
CME 



The conservation of the electromagnetic current re- 
quires the introduction of the Bardeen counterterm into 
the action. In AdS/QCD models of the CME, this typ- 
ically leads to a vanishing result for the electromagnetic 
current [30, 1321 . The problem is related to the difficulty 
of introducing a chemica,l potential conjugated to a non- 
conserved chiral charge [30, Isll • It is possible to modify 



Our main result is (|40p , which gives the chiral magnetic 
conductivity kb for an anisotropic plasma. It explicitly 
shows the dependence on the momentum anisotropy Sp. 
We also computed the CME coefficient in the holographic 
dual model and found numerical agreement with the hy- 
drodynamic result for small anisotropies. 
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APPENDIX 

Appendix A: Computation of d^u>'^ and d^B'^ 



In the following wc will use the identities 



(Al) 
(A2) 



To find an explicit expression for S^w^, we compute the 
term Lo^d^T^^'^ in two ways. First, using the hydrody- 
namic equations, we get 

u^d^.T^" = Lo,F''^'J^ = pcj.F^^M^ = pu:,E" . (A3) 

Next, substituting the stress-energy tensor ^ in this ex- 
pression, we find 

= (e + FT)w^w.a^M" + Lo^gt^^d^Pr - Atu^vt^d^v" 

= -{e + PT)u>'u''dpi^^ + uj^df.Pr - Auj^d" Inr 
- v'^uj^yi'dpA - Av'^u^df^v" . (A4) 

Using the identity (|A1I) . we find 

2 



a^c^^ = 



-u^id^PT-Adf^hiT-pE^ 



e + Pr 

~ v^v^'d.A - Av^.d.v") . (A5) 

Similar manipulations of the term B^dfj^T'^'^ lead to 

B^d^T'^'^ - B.^F'^^j^ = pBpE'^ , (A6) 

S.S^T'^'^ = -{e + PtW^u^O^B, + B>'d^PT 

- AB^v^'df.v" - B^v''v>'df,A - AB^v'^df.v'' 
= -{e + PrKdi^B^' - 2w'^£;^) - AB^a^lnr 

- B^v^yf'df.A - AB^v'^d^.v^' , (A7) 



where we used (|A2[) . From (|A6|) and (IA7|) we obtain the 
following expression: 



d,B^ = -2w^^„, - 



Bt' 



e + Pi 



id^^Pr ~ Adf^lnr - pE^ 



v^v^d.A - Av^d^v-") . (A8) 



The last two terms in (jASp and (|A8I) vanish provided 
the fluid satisfies 



a.w'^=0, w''(9„A = 0. 



(A9) 



Then (|A5|) and (jASp become identical to the expressions 
inCSl). 



Appendix B: Computation of the transport 
coefficients ^ui and ^b 

In this appendix we compute the conductivities ^^j 
and ^B by solving the system of Eqs. (f22|) - (f25l) . Fol- 
lowing [14l |. we change variables from Inr, p, T to Inr, 



p, = p/T and Pt- From p^ and ([T^ . we derive the 
thermodynamic expressions 



(Bl) 
(B2) 
(B3) 



'dp\ 


e-I-Pt 


.^^/Pt, Inr 


pT2 ' 


^aPT\ 


e + Pr 


v57^y,,in. 


T ' 


91nT\ 


1 e + Pr 


dT J,,P. 


A T 



Using 



dD dD dD 



dDf 



ODt 



dDf 



d^DB = -T^d^Pr + ^dp~p+ —^9^ In r , (B5) 
oPt op a Inr 

the first two equations, p2|) and (|23| . can be rewritten as 

dD 2D DDb Db , , 

= , TTTe -7^ = , (B7) 



OPt e + Pt 



OPt e + Pt 



'^ '^^ ^0, l^ + ^^^O. (B8) 



ainr € + Pt 



dluT e + Pj 



Note that (|B7p and (jBSP are related by the thermody- 
namic identities (jB2p and (|B3p . Using the ansatz 



D = T^d{p,\nT), Db =TdB{pMT), (B9) 

and (jBip . we obtain two differential equations from ([M)) 
and ((25|l . 

0=^^-2D. + ^ 
e + Pt T 

= T (df^diil, In r) - 2dB (m, In r))) , (BIO) 



= 



7^ -Cm 



e-hPT P 
= dfidB{pMT)-Cp. (Bll) 

These equations can be integrated to give 
1 



dBipMr) = ^Cp^ +l3{hiT) 
1 



(B12) 



d{p, In r) = ^(7/2^-^ 2/I/3(ln r) + 7(ln r) , (B13) 

where /3(lnT) and 7(lnT) are arbitrary functions of Inr. 
Substituting this back into ([24]) , p5|) , we get the conduc- 
tivities 

^. ^ C ( p^ -'^--^) +2T^Pilnr) 



ie + Pr 



2pT^ 
e + Pr 



(2^/3(lnT)+7(lnT)) 



iB=C[p~ l--^) - -^/3(lnr) . (B14) 



2 e + Pr J e + Pr 
The function 7 (Inr) is forbidden by CPT invariance [3 
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Appendix C: First-order corrected background 
geometry 



We still need to determine the constants Ci. Using 
CS]) . we replace Af in (|C3|) and obtain 



In this appendix we compute the first-order corrections 
to the background ([CT|) using the ansatz (|551) . The com- 
putation follows that for the three-charge STU model |27| 
presented in [51| and |20 |. 

We begin by substituting the ansatz ([55)1 into the equa- 
tions of motion (|44p and P5)) . We denote the resulting 
Maxwell equations, Eqs. (|45|) by A/^ (a = l,...,ri,) and 
the components of the Einstein equation, Eqn. ([^^ by 
£^A/wM,iV = 0,...,4[a;*^ = (i,.Ti,x2,a;3, r)]. Then, from 
g'^'E^tj + g'^^'Eri = 0, we find dtUi = 0, and Eu, Ert, Err, 
Ett, Mf , and M° are solved by 



The remaining equations are Eij , i^t^ , A'/° . 
From Eij we get 



(CI) 



dr r'f{r)dr 



9ij{f) 



2,r'^{diUj+djUi). (C2) 



From Eti we get 
fir) (2 , u;^(r) 



/(r) \r WT{r) 
r 2zi;i(r) 



3/M 



E<w'-6 



Va=l 



9u{r) 



5;M+3i'.(r)=4E<(r)ir(r), 



a=l 



(C3) 



where a prime denotes the partial derivative dr with re- 
spect to r. 

From Adf we get 



dr 



WL{rY'\[f{r)Af~Mr)Al- 



= a. ( -SabcAlA'^^e'^'' {djUk) + 5a6cA[;e'^"'^ [djAD 



a.Q^M 



(C4) 



Equation (|C2p depends only on gij(r) and can easily 
be solved. The integration of (jC4[) leads to 



(C5) 



WL{rY'^(rf{r)Af{r)~r~gu{r)A'^^'{r) 
^Q-{r)+Ct. 



Here C° are some integration constants, which can be 
fixed as 

Cf - -g?(r+) - C.^L(r+)i/V+Ag'(r+) , (C6) 

with r+ as in (|59p and C^ = gti{r^). This can be solved 
for ilf(r), 



i^M 



dr - 



rj[r')WL(r)' L 
Qr+^S'(r+)u;i(r+)V2 + r'g,,(r')Ag'(/)' 



(C7) 



/'(r) /2 , ii;^(r) 



/(r) Vr WT{r)J 3/(r) 



2wL(r) 



5^.(0 +5*. W 



1 



v,a=l 



WL{rY/'^rf[r) 



9u{r) 



(C8) 



where 



/(r) = E4<W(Q?W-Q?(^+) 



a=l 



ar+t.i(r+)i/2^g'(r+) 



(C9) 



A homogeneous solution of this equation gu(r) = 
g).^ (r) = f{r) can be generated by the infinitesimal co- 
ordinate transformation 

dr 



dt ^f dt — e{dx + dy + dz) , dz -^ dz + e 

dv dv 

dx -^ dz + e^^ , dy -^ dy + e 



r'^WL 



r^WT 



r^WT 



(CIO) 



Then, using this homogeneous solution and Appendix [Dl 

[P{r) = f{r) and E{r) = rwL{ry^^ there], we bring 
to the integrablc form 



dr ( WL{ry/^rfir)d, 



fir) 



I{r). (CU) 



Solving this equation for gu (r) and fixing the integration 
constants at r+, we get 



~9u{r) = fir) / dr' 



1 



dr I{r ) 



WL(r')i/2r'(/(r'))' 
-WLir+)'/'r+f'ir+)a). (C12) 

In the Landau frame we require u^t^^'^ = 0, which 
in particular implies the absence of corrections to T**. 
Holographic renormalization |52l | translates this into a 
constraint for the r~^ coefficient of guir) which is pro- 
portional to the first correction of T*', 



lim r'^guir) =0. 

r— >-oo 



(C13) 



In the limit r — > oo, we have the asymptotics 
/(r)=0(r2), WLir)^0{\), 

dr'l{r') = 0{l), 



(CM) 

and, from the vanishing of the r~^-coefficient of guir), 
we obtain the following equation for Ci: 

/•OO 

WLir+)'/\+f'ir+)a = / dr'l{r')=I, + I2 • C, , 

(C15) 



where we defined the integrals 

11^4/ dr'^<(r')(Q,"(r')-g?(r+)) 

Jta- 1 



"'+ 0=1 



+ 2 SabcA'^^{r+)Al{r+)e^=^ {d,Al) (C16) 



-Sa,cAl{r+)Al{r+)Al{r+y^'^{d,Uk) 



and 



I2 ^ 4 / dr' 5] Ag'(r') (-«;L(r+)i/V+yi;5'(r+) 

"''■+ a=l 

= 4^,.(r+)i/2^+ ^^f5(,^)^-(^^) . (C17) 



Solving this for Ci, we eventually get 
4 






+ -Sa,cAl{r+)Al{r+)e^=^ {d,Al) 



(C18) 



Appendix D: Integrable form of a linear ordinary 
differential equation 

In this appendix we present a method to bring an arbi- 
trary linear ODE of second order to an integrable form. 
Let us consider a general form of this equation 



G{g", g\ g, r) = g"{r) + a{r)g'{r) + b{r)g{r) = c{r). 



(Dl) 
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If we know a homogeneous solution P{r) of this equation, 
i.e. 

G{P",P',P,r)^0, (D2) 

then we can make the substitution 

g{r) -> P{r)Q{r), Q' {r) -> u{r) (D3) 

and lower the order of the differential operator (|D1I) 

P'{r)- 



G = P{r) u'{r) + 



a(r)+2. 



P{t) 



u{r) 



= P{r){u'{r) + F{r)u{r)). (D4) 

The term in the brackets can be represented as 

1 



u'{r) + F{r)u{r) 



A{r) 



dr {A{r)u{r)) , (D5) 



where 



A{r) = exp <^ / F{r)dr \ = P{rY cxp <^ / a{r)d' 



(D6) 

Taking into account (jD3p . we finally bring ()D1[) to the 
following integrable form 



1 



dr P{rYE{r)dr 



9ir) 



P{r)E{r) V V^(^) 

where we defined 

E{r) = exp < / a{r)dr 



c(t) . (D7) 



(D8) 
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